1 

UTHEP-342 
UTCCP-P-15 
August 1996 

Phase structure of QCD for general number of flavors * 

Y. Iwasaki?< b K. Kanaya a < b S. Kaya b S. Sakai^ and T. Yoshie a < b 

a Center for Computational Physics, University of Tsukuba, Ibaraki 305, Japan 

b Institute of Physics, University of Tsukuba, Ibaraki 305, Japan 

c Faculty of Education, Yamagata University, Yamagata 990, Japan 

We investigate and elucidate the phase structure of QCD for general number of flavors Nf with Wilson quarks, 
varying Nf from 2 up to 300. Based on numerical results combined with the result of the perturbation theory we 
propose the following picture: When Nf > 17, there is only a trivial fixed point and therefore the theory in the 
continuum limit is trivial. On the other hand, when 16 > Nf > 7, there is a non-trivial fixed point and therefore 
the theory is non-trivial with anomalous dimensions, however, without quark confinement. Theories which satisfy 
both quark confinement and spontaneous chiral symmetry breaking in the continuum limit exist only for Nf < 6. 
We also discuss the structure of the deconfining phase at finite temperatures for the small number of flavors such 
as Nf — 2 and 3, through a systematic study of it for general number of flavors. 



1. Introduction 

In our previous work [Q, it was shown that in 
the strong coupling limit, when the number of 
flavors Np is greater than or equal to 7, quarks 
are deconfined and chiral symmetry is restored at 
zero temperature, if the quark mass is lighter than 
a critical value. On the other hand, it is well- 
known that when Np exceeds 16^ , asymptotic 
freedom is lost. 

In this work we investigate the problem of what 
is the condition on the number of flavors for quark 
confinement and spontaneous chiral symmetry 
breaking in the continuum limit, using the Wil- 
son quark action and the standard one-plaquette 
gauge action. We first clarify the phase structure 
at zero temperature for general number of flavors, 
in particular, for Np > 7. When the phase dia- 
gram becomes clear, we are able to see what kind 
of continuum limit exists and eventually answer 
the problem given above. 

In order to reveal the structure of the deconfin- 
ing phase, we increase Nf up to 300, because the 
region of the deconfining phase at zero tempera- 
ture becomes larger with larger Np, and conse- 
quently the structure becomes clear, as shown be- 
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low. Then combining all data including the case 
of small Np, we are able to conjecture the phase 
structure for general number of flavors. Through 
this systematic study of the deconfining phase for 
general number of flavors, we are also able to un- 
derstand the behavior of the quark mass m q (de- 
fined through an axial-vector Ward identity) and 
m 2 as functions of 1/K in the deconfining phase 
at finite temperatures in particular at small (3 for 
small number of flavors such as Np — 2 and 3. 

2. Simulation parameters 

The lattice sizes are 8 2 x 10 x N t (N t = 4, 6 or 
8), 16 2 x 24 x N t (N t = 16) and 18 2 x 24 x N t 
(Nt — 18). We use an anti-periodic boundary 
condition for quarks in the t direction and peri- 
odic boundary conditions otherwise. When the 
hadron spectrum is calculated, the lattice is du- 
plicated in the direction of lattice size 10 for 
N < 8, which we call the z direction. We call 
the pion screening mass simply the pion mass and 
similarly for the quark mass. The current quark 
mass rriq is defined by an axial-vector Ward iden- 
tity §,§. 

We use the hybrid R algorithm for the quark 
matrix inversion. As Np increases we have to 
decrease At, such as At =0.0025 for Np = 240, 
to reduce 0(At 2 ) errors. We have checked that 
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Figure 1. The effective pion mass obtained at z=7 
and the current quark mass for a free quark versus 
1/K for Nt — 4 and 8 with periodic and anti- 
periodic boundary conditions on 40 3 x Nt lattices. 

the errors are sufficiently small selecting typical 
cases. 

3. Free case: (3 = oo 

We first investigate the free case, because the 
free case is important as a reference to understand 
the structure of the deconfining phase. We will 
show later that the 1 /K dependence of the pion 
mass in the deconfining phase is quite similar to 
that of free quarks for any Np and (3. 

Because the Wilson formalism lifts the doublers 
by the Wilson term, there is only one pole in the 
free quark propagator at 1/K ~ 8. However, 
there is reminiscence of doublers. 

In Fig. [| we plot the effective pion mass ob- 
tained at z — 7 from the pion propagator. The 
four cases shown there correspond to those for 



Nt — 4 and 8 with an anti-periodic or the peri- 
odic boundary condition in the t direction. 

In the case of the periodic boundary condition, 
there are almost no differences between the Nt — 
4 and 8 cases. The vanishing of m q at 1/K = 8 
and the behavior of for 1/K > 8 roughly given 
by ~ — 8) are all as naively expected. 

For 1/K < 8, the behavior of m% is somewhat 
complicated. The vanishing of at 1/K = 4 
is due to doublers; one of momenta in spatial di- 
rection is 7r. At 1/K ~ 6, the momentum-energy 
(dispersion) relation is completely opposite to the 
usual one: The energy is a decreasing function of 
the magnitude of momentum for the case of equal 
momenta in three spatial dimensions. The energy 
becomes infinity at zero momenta like the infinite 
mass case, while it is finite for other modes of mo- 
menta. Therefore the appearance of the peak of 
run at 1/K ~ 6 is due to a superposition of con- 
tributions from several modes of momenta. The 
behavior of m q defined through the axial-vector 
Ward identity is also complicated at 1/K < 8 
due to doublers. In the following, we refer to the 
rightmost m q = state (at 1/K — 8 for j3 = oo) 
as simply the massless quark, or m q = 0, dis- 
tinguished from the other m q = states, unless 
otherwise stated. 

When an anti-periodic boundary condition is 
imposed, the position as well as the height of 
the peak of at 1/K ~ 6 does not change so 
much. On the other hand, for the massless quark 
at 1/K — 8, the pion mass takes a value which 
is roughly twice the lowest Matsubara frequency, 
for the both cases of Nt = 4 and 8. The behavior 
around 1/K — 4 is similar, but indicates more 
complicated structure. 

4. Kd versus Np 

We now study the case of j3 — as in the pre- 
vious work [|J. Fig. ^a shows the results of the 
deconfining transition points Kd at [3 = ob- 
tained on Nt = 4 lattices for various numbers of 
N F (see also Fig. §a). 

The data of plaquette (Fig. |) for N F = 240 
at /3 = indicates that the deconfining transi- 
tion is of first order and that the location of the 
transition point is independent of Nt for large Nt 
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Figure 2. a) The transition point 1/Kd at j3 = 
versus ATp at iVj = 4 and N t > 8. For clarity, 
data at JV t = 8 for N F = 300 is slightly shifted 
in the figure. Shaded region was investigated in 
our previous study M. b) Plaquette versus l/K 
at (3 = for various number of -ZVj? at iV t = 4. 

(1/ATd ~ 8.1 for N t = 4 and 1/A d ~ 7.75 for 
N t — 8 and 16). Therefore we conclude that the 
transition is bulk as is confirmed in our previous 
work for Np = 7 jjj. In Fig. ||a the transition 
points at N t > 8 for Np = 7, 240 and 300 are 
also included. These values are roughly those for 
the bulk transition points at zero temperature. 

As Np increases up to 240, 1/Kd approaches 
toward or exceeds the value 8, which is the value 
of 1/K for a massless free quark. Because of this, 
for Np > 240, we are able to see a wide region 
of the deconfining phase in the whole range of /?. 
Therefore we first intensively investigate the case 
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Figure 3. Plaquette at = for Np = 240 at 
N t = 4, 8, and 16. 

Np = 240, and then decrease Np. 

5. Np = 240 

Fig. [| shows the results of and 2m q for 
Np = 240 at = 0, 2.0, 4.5, 6.0, and 100 on 
the Nt = 4 lattice. A very striking fact is that 
the shape of and 2m q as a function of 1/A 
only slightly changes for 1/K <8 when the value 
of P decreases from oo down 0. Only the position 
of the local minimum of at 1/K ~ 8, which 
corresponds to the vanishing point of m q , slightly 
shifts toward smaller 1/A. We obtain similar re- 
sults also for Nt — 8. 

From these data we propose the phase diagram 
in Fig. H for Nf = 240. The dark shaded line is 
the phase boundary between the confining phase 
and the deconfining phase at zero temperature. 
When Nt = 4 or 8, the boundary line bends 
down at finite (3 as shown in Fig. ||, due to the fi- 
nite temperature phase transition of the confining 
phase. 

The dashed line corresponds to the m q = 
line. This line also corresponds to the minimum 
point of mj. We have also calculated the quark 
propagator in the Landau gauge and checked that 
chiral symmetry of the propagator, 75(7(2)75 = 
—G(z), is actually satisfied on the m q = line. 

The results for the case Np = 300 are essen- 
tially the same as those for Np — 240 except for 
very small shifts of the transition point and the 
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Figure 4. and 2m q versus 1/K for Nf = 240 
at iV t = 4. 

minimum point of m^. 

6. Renormalization group flow 

The m q = point at (3 = oo is the trivial IR 
fixed point for Np > 17. The phase diagram 
shown in Fig. || suggests that there are no other 
fixed points on the m q = line at finite (3. 

In order to confirm this, we investigate the 
direction of Renormalization Group (RG) flow 
along the m q — line for Np = 240, using 
a Monte Carlo Renormalization Group (MCRG) 
method. 

One problem here is that it is practically im- 
possible to make simulations in the massless limit 
at zero temperature due to the existence of zero 
modes in the quark matrix. Therefore, we im- 
pose an anti-periodic boundary condition in the t 
direction. 

We make a block transformation for a change 
of scale factor 2, and estimate the quantity A/3 
for the change of a — > a' — 2a. We generate 
configurations on an 8 4 lattice on the m q = 
points at (3 — and 6.0 and make twice blockings. 
We also generate configurations on a 4 4 lattice 
and make once a blocking. Then we calculate A/3 
by matching the value of the plaquette at each 
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Figure 5. Phase diagram for Np — 240. 
step.Q 

From the matching, we obtain A/3 ~ 6.5 for 
(3 — and 10.5 for (3 = 6.0. The value obtained 
from the perturbation theory is A/3 ~ 8.8 at j3 = 
6.0 for Np — 240. The signs are the same and the 
magnitudes are comparable. This suggests that 
the directions of RG flow on the m q = line at 
(3 = and 6.0 are the same as that at (3 = oo for 
Nf — 240. This further suggests that there are 
no fixed points at finite (3. All of the above imply 
that the theory is trivial in the case of Np = 240. 

7. 240 > N F > 17 

Now we decrease the value of Np from 240. 
Fig. H shows to \ and 2m q for various Np, from 
300 down to 18, at (3 = 6.0 and 0. 

When (3 — 6.0, the shapes of are almost 
identical to each other, except for a slight shift 
toward smaller 1/K as Nf decreases. 

When (3 — 0, the boundary of the first order 
phase transition between the deconfining phase 
and the confining phase moves toward smaller 
1/K and therefore the range of the deconfining 
phase decreases: I/K4 ~ 8.5, 8, 7.6, 7.2, 6.1 for 
N F = 300, 240, 180, 120 and 60, respectively, as 
shown in Fig. ^. For Np — 18, the value of 
decreases down to 5.0. 

Due to the fact that the confining phase invades 

1 It is known for the pure SU(3) gauge theory, in par- 
ticular in the deconfining phase, that one has to make a 
more careful analysis using several types of Wilson loop 
to extract a precise value of A/3. We reserve elaboration 
of this point and a fine tuning oi 1/K at each fi for future 
works. 
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Figure 6. ml and 2m q for N F < 18 ,N t = 4 at 
/3 =6.0 and 0.0 versus 1/K. 

the deconfining phase, the massless line in the 
deconfining phase hits the boundary at finite (3 
when Np becomes small. For example, in the 
case of Np — 18, as shown in Fig. [j], it hits at 
(3 = 4.0 — 4.5. 

Although the area of the deconfining phase de- 
creases with decreasing Np, the shape and the 
position of m\ in the part of deconfining phase 
only slightly change from Np = 300 to 18. The 
values of m q as functions of 1/K also show only 
slight changes in the deconfining phase. These 
facts, combined with the perturbative result that 
for Np > 17, (3 — oo is the IR fixed point, sug- 
gest that the RG flow along the massless quark 
in the deconfining phase is the same as that for 
Np — 240. In this case, the theory is trivial for 
N F > 170 

2 Although we have investigated only down to the Np = 
18 case, we do not expect any qualitative differences be- 
tween the cases of Np = 17 and 18. 



0.2 

K 

0.1 



In=18 


deconfining phase _ 








■ confining phase 




; -m (N=4) = 
o K(N=4) 

■ - -. - m (N =8) = 
□ K(N=8) 


K 

"~d N t large 
_i — , — , — — , — , — , — i — , — , — , — 



p 

Figure 7. Phase diagram for Np = 18. 



8. 16 > N F > 7 

When Np < 16, the theory is asymptotic free. 
On the other hand, quark confinement is lost for 
Np > 7 even in the strong coupling limit /3 = 
jlj]. Therefore the question is what happens for 
the cases 16 > Nf > 7 in the deconfining phase. 

We have intensively simulated the cases N F = 
12 and 7. See the results for the case Nf = 7 in 
Fig. ||. The gross feature of the phase diagram 
does not much differ from the case of Np = 18 
which is shown in Fig. [?]. 

One crucial point is that there is no chiral limit 
in the confining phase at least for (3 < 4.5. This 
suggests that quarks are not confined in the con- 
tinuum limit. 

A more subtle point is whether there is a dif- 
ference between the case Np > 17 and the case 
16 > Np > 7 concerning the deconfining phase. 
For Np < 16, the massless point at (3 — oo is 
a UV fixed point. Therefore the direction of the 
RG flow around (3 = oo is opposite to that for 
Np > 17. In order to study the phase struc- 
ture around j3 = 0, we have performed a MCRG 
study for the case Np = 12 similar to that for 
the case of Np — 240. We have investigated the 
RG flow along the maximum of m\ , because this 
line roughly corresponds to an infinite quark mass 
line in a sense explained in Sec. |[ We have found 
that the RG flow at f) = and (3 = 6.0 are op- 
posite and the flow changes the direction around 
(3 — 4.0. This result suggests that there is a fixed 
point on the m q = line. The flow of RG on 
the m q = line at (3 = 5.0 is the same as that 
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at = oo, which is opposite to that of the case 
of Np — 240. This implies that the location of 
the fixed point is at (3 < 5.0. Because of the fact 
that the massless line m q = hits the bound- 
ary between the two phases around (3 = 4.0, it 
is difficult to determine the location of the fixed 
point. Although we have not been able to lo- 
cate the fixed point, we propose the phase struc- 
ture, based on the consideration given above, that 
there is a non-trivial IR fixed point at finite (3 Q 
for 16 > Np > 7. Of course we have to make 
a detailed RG study to confirm that there is a 
non-trivial fixed point, which we reserve for fu- 
ture works. The existence of a non-trivial fixed 
point implies that the theory in the continuum 
limit is a non-trivial theory with anomalous di- 
mensions, however, without confinement. 

9. N F < 6 

For Nf = 2 and 3 at finite temperatures, it was 
shown that m q has an unexpected 1/K depen- 
dence in the deconfining phase when (3 is smaller 
than about 5.3 The as a function of 

1/K also shows a cusp. In this section we concen- 
trate on these issues. (For the discussion on the 
phase structure at zero temperature, see Ref. || 
and references in it.) 

We point out that the behaviors of m q and 
TOtt in the high temperature deconfining phase 
for Np < 6 at (3 < 5.3 are similar to those in the 
deconfining phase at zero temperature for larger 
Np. In Fig. ||, the value of 2m q is plotted vs. 1/K 
for from N F = 240 down to N F = 2 at ft = 4.5. 
Results at (3 — 5.0 are similar. The value of 2m q 
for the free quark case is also plotted. There are 
two major trends when Np decreases from 240 
down to 2: One is an overall shift of the function 
2m q in the deconfining phase to smaller 1/K and 
shrinkage of the shape with decreasing Np. The 
other is a shift of the transition point from the 
confining phase to the deconfining phase toward 
smaller 1/K. 

It should be noted that except for the overall 
shift and shrinkage mentioned above, the shape of 
the function 2m q in terms of 1/K in the decon- 
fining phase is quite similar for all Np and also 
similar to that of a free quark. Thus the behavior 



of m q as well as m^. as functions of 1/K in the 
deconfining phase at /?<5.3 can be understood 
essentially as those of free quarks for 1/ K < 8. 

We further note that for Np — 2 and 3, the 
massless quark point (the rightmost m q — 0) in 
the deconfining phase is different from the K c 
point in the confining phase for /? < 5.3. The 
difference is illustrated in the phase diagram for 
N F = 2 on the (0, K) plane (see Fig. 0). On 
the K c line, both m^. and m q vanish in the con- 
fining phase at zero temperature The finite 
temperature transition for Nt = 4 occur on the 
line K t . The massless quark line for Nt — 4 in the 
high temperature deconfining phase agrees with 
the K c line at f3 ~ oo. However, they deviate 
from each other at (3 < 5.3. This corresponds to 
the fact that the value of m q does not depend on 
the phase when (3 > 5.5, while it strongly depends 
on the phase when /3<5.3 ]§J|. 

10. Conclusions 

Based on numerical results combined with the 
result of the perturbation theory we propose the 
following picture: There are three categories de- 
pending on the number of flavors Np: a free the- 
ory for Np > 17, a non-trivial IR fixed point for 
16 > Np > 7 and confinement and spontaneous 
chiral symmetry breaking for Np < 6. The con- 
clusion that there are three categories depending 
on Np is in accord with recent theoretical works 
of N = 1 super-symmetric QCD |Tc[ |. 

We also find that, for any Np > 7 and for Np < 
6 at (3 < 5.3, the 1/K dependences of m q and m 2 
in the deconfining phase are essentially the same 
as those of free quarks for 1/K < 8, which show 
non-trivial structure due to doublers. The pecu- 
liar behavior of m q as well as the cusp structure of 
in the deconfining phase at f3 < 5.3 for Np = 2 
[|| can be understood in terms of the overall 
shift in 1/K and shrinkage of this structure. 

Numerical simulations are performed with HI- 
TAC S820/80 at KEK, and Fujitsu VPP500/30 
and QCDPAX at the University of Tsukuba. This 
work is in part supported by the Grants-in-Aid 
of Ministry of Education, Science and Culture 
(Nos.07NP0401, 07640375 and 07640376). 
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Figure 8. m\ and 2m q for Nf = 7 versus 
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Figure 10. Phase diagram for Np = 2. 



